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The present invention relates to a method of measuring the 
dilution of phase modulate d spins as well as to a magnetic rep o- 
naiice^ maging device which can be used for performing this 
method^ 

The present invention finds particular application in regard 
to medical diagnostic magnetic resonance imaging, but it is to 
be appreciated that the present invention also finds application 
in magnetic resonance spectroscopy and magnetic resonance imag- 
ing for other applications. 

In biological tissues water protons are often compartmental- 
ized in terms of their molecular mobility. Commonly, the two- 
pool model is used to explain macroscopic relaxation properties 
resulting from this compartmentalization : One pool, usually the 
largest in size, is associated to the "free" mobile bulk water, 
whereas a second pool is associated to motional restricted pro- 
tons bound to macromolecules such as proteins or lipids. In 
"Magnetic Resonance in Medicine" (Vol. 10, P. 135-44 (1989)) S. 
D. Wolff and R. S. Balaban first coined the terms free and re- 
stricted proton pools to describe these different pools. Most 
importantly, both pools are coupled by an exchange of magnetiza- 
tion via chemical exchange and dipolar coupling. This phenome- 
non, commonly termed magnetization transfer (MT) , is described 
in U.S. Pat. No. 5,050,609 and how to make use of it in magnetic 
resonance imaging is reviewed in "Magnetic Resonance Quarterly" 
(Vol. 8, P. 116-37 (1992) ) . 

In medical diagnostic imaging the quantity of the restricted 
(macromolecular) proton pool is highly desirable as it is ex- 
pected that this pool directly represents tissue structure and 
tissue integrity. However, depicting of the restricted proton 
pool is not possible because in most tissues the transverse mag- 
netization of these protons decays with a time constant of ap- 
proximately 10 fxs or less. Unfortunately, conventional magnetic 
resonance imaging systems are not capable of sampling such ul- 
trashort signals. Therefore, the restricted proton pool can be 
depicted only indirectly by exploiting the MT phenomenon. 

Current methods for determining relaxation parameters of the 
two-pool model including the relative proton density employ both 
the formalism of the coupled Bloch Equations and spectral selec- 
tive RF pulses with the intention to irradiate either the re- 
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stricted or the free pool. Following spectral selective radio 
frequency (RF) irradiation, the signal response of the free pool 
is sampled and related to the model parameters . 

In "Magnetic Resonance in Medicine" (Vol. 29, P. 759-766 
(1993)) Henkelman R. M. et al . have derived an equation which 
relates the attenuation of the steady state magnetization to the 
resonance offset and power of a continuous wave saturation pulse 
and to some fundamental model parameters including the relative 
size of the restricted pool. The pool parameters can be obtained 
then by fitting this equation to several measurements obtained 
with different RF powers and frequency offsets. So far, this 
method could not be applied in a clinical setting because it ex- 
ceeds the current limit by the specific absorption rate (SAR) 
and because measurement time is impractically long. Addition- 
ally, this technique suffers from the constraint of an addi- 
tional measurement of the apparent relaxation time Tl . Moreover, 
prior knowledge regarding the lineshape and the transverse re- 
laxation time of the restricted pool is required. 

Instead of off -resonant saturation pulses, another method, 
presented by Daniel Gochberg et al . in "Magnetic Resonance in 
Medicine" (Vol . 41, P. 1065-1072 (1999)), employs a train of 
short on-resonant inversion pulses. This approach benefits from 
a low SAR as the inversion pulses are separated by 120 ms. How- 
ever, it suffers from a long measurement time as well. If more 
than one slice is required to be imaged, measurement time scales 
linearly with the number of slices as this method is not capable 
of true multislice imaging. So far, no in vivo results have been 
generated with this method. 

More recent work, presented by John Sled et al . in "Magnetic 
Resonance in Medicine" (Vol . 46, P. 923-931 (2001)) and by Vasily 
Yarnykh in "Magnetic Resonance in Medicine" (Vol . 47, P. 929-939 
(2002)) is also based on the steady state solution of the cou- 
pled Bloch Equations but regards pulsed RF saturation by incor- 
porating numerical calculations. However, these newer methods 
are still not suited for routine use due to long acquisition 
times and extensive calculations. 

It is an object of the present invention to provide a above 
described method which is capable of eliminating the drawbacks 
of the prior art and which enables a fast and low specific ab- 
sorption rate (SAR) method for measuring the macromolecular pro- 
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ton density. 

The present invention relates to magnetic resonance spec- 
troscopy or imaging techniques in which an initial preparation 
of the proton magnetization is performed by two radio frequency 
pulses which causes all, or a portion of, the proton magnetiza- 
tion in the volume of interest to be oriented longitudinally. 
The application of a gradient field between these RF pulses 
causes the longitudinal magnetization to be modulated along the 
field direction. Such a preparation scheme, commonly known as 
stimulated echo preparation, does only affect protons of the 
free mobile tissue water. The restricted proton pool is not af- 
fected, because a RF pulse separation time is chosen that is 
much longer than the transverse relaxation time of the re- 
stricted proton pool. The present invention, therefore, employs 
above preparation scheme such as to label the free water protons 
and uses these protons as an intrinsic indicator to measure the 
size of the restricted pool. Following labelling, the concentra- 
tion and therefore the signal intensity of the labelled protons 
will decrease as the labelled protons "dilute" into the re- 
stricted proton pool via magnetization transfer. Additionally, 
the labelled protons are also subject to longitudinal relaxa- 
tion. Probing the labelled magnetization by means of a third RF 
pulse after different mixing times yields an indicator decay 
curve. From this curve the relaxation and dilution effects are 
separated by a bi- exponential analysis and the size of the re- 
stricted proton pool is calculated according to indicator dilu- 
tion theory. 

In accordance with a preferred embodiment of the invention, 
an additional RF pulse is applied between the labelling pulses 
and the final readout pulse. This additional pulse is intended 
to change the magnetization in the free proton pool only. Addi- 
tionally, a mixing time, i.e. the time between the second and 
the last RF pulse, is used that is twice the time needed for 
both proton pools to restore equilibrium. Then, the size of the 
restricted pool is determined from two acquisitions, where the 
flip angle of said RF pulse is 0° in one run and 180° in the 
other run. Said RF pulse can also be a composite pulse that, pro- 
vides a constant specific absorption rate independent of the ef- 
fective flip angle. 

Alternatively, in a variation of the method, the flip angle 
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other suitably equipped MR I system that may be programmed to 
measure the macromolecular proton concentration. As illustrated, 
a magnet 10 generates a static, fundamental magnetic field along 
a z axis 12, in which an object or the body of a patient 14 to 
be examined is situated. The system additionally comprises gra- 
dient amplifiers 16, gradient coils 18 , transmitter 20, RF 
power amplifier 22, RF coils 24 for generating pulse sequences 
for application to selected slices of the patient's 14 anatomy 
or the sample. The control of the pulse sequence is done by a 
sequence control unit 26, which can be programmed through a scan 
control interface 28. Since software techniques for generating 
pulse sequences with the characteristics defined below are be- 
lieved to be well-known to those skilled in the art, such pulse 
generating techniques will not be described in further details 
herein. The signal generated by the pulse sequence is received 
by signal receiver 3 0 and digitized at digitizer 32 for applica- 
tion to an arithmetic unit 34 for processing in accordance with 
the technique of the invention. The processed signal is then 
displayed on a display unit 36. Data storage 38 and filming with 
a camera device 40 may be provided additionally. To synchronize 
the pulse sequence with physiological signals from patient 14, a 
cardiac synchronization unit 42 and a controlling device 44 for 
respiratory motion may also be provided. 

Now, suppose water protons in tissue to be in two different 
states in regard to their molecular mobility. The compartmen- 
talization in terms of molecular mobility is illustrated in Fig. 
2. This figure shows a simplified model of relaxation in het- 
erogeneous tissues. This is commonly know as the two-pool model. 
The pool A corresponds to X H spins in "free" mobile bulk water, 
whereas the second pool B, usually much smaller in size, corre- 
sponds to motional restricted X H spins bound to macromolecules . 
Pool B will be referred to as restricted or macromolecular pool. 
Each pool is characterized by its intrinsic relaxation rates R 2 
and R 2 and by its size M a ° and M b ° . Most importantly, there is an 
intermediate to fast exchange of magnetization between both 
pools as indicated by the first order forward and backwards 
transfer rates k f and k b . The fundamental parameter which can be 
measured with the present invention is the size of the macromo- 
lecular pool. For simplicity the molar fraction f will be used 
in further calculations, which is the relative size of the mac- 
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romolecular pool and which can be written as 



f M » 



Since the condition for microscopic reversibility has to be ful- 
filled, M a °'k f = M b °'k b and Equation [1] can also be written as 



As known by those skilled in the art of biomedical engineering, 
a common concept to determine the unknown distribution of vol- 
umes or fractional sizes is to measure the dilution of an indi- 
cator. The indicator can be of any type, but must be inert and 
directly accessible to a measurement such as in pool A in Fig. 
3(a) for instance. After application of the indicator into pool 
A, indicated by the filled circles, an exchange process with 
pool B will start. As soon as the dilution of the indicator, 
caused by diffusion or other exchange processes, has reached a 
steady state, i.e. the concentration of the indicator is the 
same in both pools (Fig. 3(b)), the relative volume fraction or 
pool size can be calculated from the change in indicator concen- 
tration because 

[Css]=[c 0 ]j^ [3] 

where [c 0 ] is the initial concentration in pool A and [c ss ] is the 
steady state concentration. 

The present invention uses labelled spin magnetization as an 
inherent and exchangeable indicator. Spin labelling is achieved 
with a stimulated echo preparation scheme, consisting of two 
successive RF pulses, preferably with a flip angle of 90°, and a 
gradient field between. For those skilled in the art it is well 
known, that the gradient field causes a phase modulation of the 
transverse magnetization along the field direction. After apply- 
ing the second RF pulse, the phase modulation turns into an 
modulation of the longitudinal magnetization. The present inven- 
tion uses a RF pulse separation which is much longer than the 
transverse relaxation time of the restricted proton pool. This 
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ensures, that only spins in the free proton pool are labeled 
even if the RF pulses partly saturate the restricted pool. Fol- 
lowing the generation of this "indicator" , subsequent decay of 
the labelled spins with increasing mixing time is caused primar- 
ily by two processes. These are the dilution effect and Tl re- 
laxation. This process can be modelled by regarding relaxation 

and magnetization exchange in the two pool system: 

dM a {t) 



dt 

and 



-=-M a {t){R la +k f )+M b {t)k b [4] 



dM.lt) 

—^=-M b (t)(R lb ^k b )^M a (t)k f [5 ] 

where the subscript a and Jb refer to the free and restricted 
pool and M(t) is the magnetization of the labelled X H spins. 

Solving Equations [4] and [5] with the condition of M b (t=0) 
= 0 gives a biexponential function for the decay of the labelled 
spins. Only the solution for the free pool is considered here 
because the bound pool will not contribute to the measured sig- 
nal due to its extremely high R 2 : 

M a (t) = M 0 (C 1 exp(~A 1 t)+C 2 exp(-A 2 t)) [6 ] 

with 



*i, 2 = 0.5(k b +k f + R la + R lb ^ 



Ai-A 2 



[8] 



and 



£ _ \\ — k f —R la 
A-A 2 



[9] 



where M 0 is the maximum magnetization available immediately after 
labelling. The rate X x is a "fast" rate which is responsible for 
the quick approach to a steady state between both pools. In con- 
trast, X 2 is a "slow" rate and is roughly comparable to R 1 ob- 
tained from a conventional T 2 measurement. It is notable that X x 
and X 2 are identical to the rates given by the general solution 
of the coupled Bloch equations after disturbance of the equilib- 
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rium (Journal of Magnetic Resonance, Vol. 31, P. 207-229 
(1978)). However, the constants C, and C 2 found here differ from 
this solution. 

in most tissues the backward transfer rate k b is much higher 
than the longitudinal relaxation rates. At 1.5 T k b is approxi- 
mately 16 times higher in white matter, 35 times higher in gray 
matter and 70 times higher in muscle than any of the longitudi- 
nal relaxation rates (Magnetic Resonance in Medicine, Vol. 33, 
P. 476-482 (1995), Magnetic Resonance in Medicine 35:277 
(1996)). Therefore, the condition k b »R la ,R lb allows to reduce the 
constants C 1 and C 2 to: 



C '=7TT 1101 



and 



c 2 = Wl tin 

Rewriting Eq. [6] for the measured signal intensity obtained 
from a stimulated echo will result in : 

S(,) = S o -JL(/exp(-A,0+exp(-A 2 0) [12] 
where S 0 is the maximal possible signal intensity that would be 
expected for a mixing time of zero. Thus, f can be calculated 
from a bi -exponential fit to a stimulated echo data set obtained 
with several different mixing times. It is important to note 
that f then is independent of the initial state of both pools, 
ineffective labelling of the free pool will result in a lower 
dynamic range of the decay curve, but this will be reflected 
only in a decrease of S 0 and in the goodness of the fit. The re- 
sult is also insensitive to any incidental irradiation of the 
restricted pool as long as this does not occur during the mixing 
period. For a mixing period much longer than l/X, (at 1.5 T this 
corresponds to 150 - 200 ms in brain white matter) , Equation 
[12] can be rewritten as 

S(t)=S 0 ^txp(-\ 2 t) [13] . 

When compared to Equation [3] it is obvious that labelled spins 
can be treated as an indicator provided longitudinal relaxation 
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Now, consider a two-point pulse sequence such as in FIG. 4, 
which makes use of the dilution effect of phase modulated spins. 
Fig. 4 shows a preferred MR pulse sequence for the determination 
of the mac r onto lecular proton fraction. This sequence exemplarily 
employs gradients for slice selection (S) , phase encoding (P) , 
frequency encoding (R) , and labelling (M) . The sequence is per- 
formed twice without and with a 180° RF pulse placed in the cen- 
ter of the mixing period. The first run serves as reference scan 
to asses the magnitude of the relaxation term exp ( -A 2 TM) which is 
equal for both runs. The intention of the second run is to in- 
troduce an imbalance in the system. This is achieved with the 
180° inversion pulse which changes the sign of the phase of the 
spins in the free pool only. This will be followed by a dilution 
of spins with opposite phase to the restricted pool and by a di- 
lution of originally labelled spins from the restricted pool to 
the free pool. Immediately before the 180° RF pulse the signal 
from the labelled spins is 

S'=S 0 jr^GQ(-A 2 7M x ) [14] . 

The signal intensity achieved in the first run is 

S,=S'exp(-A 2 73lf 2 ) [is] 

and for the second run with the 180° RF pulse the signal inten- 
sity is 

S 2 =-S'y^exp(-A 2 rM 2 )+^ [16] 

where the second term in Equation [16] gives the dilution from 
originally labelled spins in the restricted pool to the free 
pool. In above equations the inversion pulse is assumed to be 
perfect. However, the inversion pulse may also affect the spins 
in the restricted pool which is considered by the parameter p 
that gives the relative saturation of the restricted pool; total 
saturation will result in a p value of 0 while no saturation 
will be indicated by a value of 1. Usually, the restricted pool 
is expected to be smaller then the free pool (f<l) . Then, if the 
signal intensity is obtained from a magnitude image the molar 
fraction is given by 

, _ (s-s 2 ) 

J ~(pSi+S 2 ) Cl7] - 
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If R 2b is known or can be estimated, p can be obtained from nu- 
merical simulations of the effect of the inversion pulse on the 
bound pool. Equation [17] can be implemented to calculate pixel - 
by-pixel parameter images showing the molar fraction. To measure 
the absolute macromolecular proton concentration, the proton 
concentration of the free water pool has to be determined. This 
can be achieved, for instance, by estimating the size of the 
long T 2 -component from an additional multi-echo experiment and by 
scaling it to a reference water sample of known temperature. 

The pulse sequence of Fig. 4 was implemented on a Intera™ 
1.5T whole-body scanner. 

The pulse sequence and method was validated with phantoms 
containing known concentrations of agar gel and bovine serum al- 
bumin. Additionally, the pulse sequence was used to evaluate the 
macromolecular content in brain tissue of several volunteers. 
Fig. 5 illustrates an application of the method of mapping the 
relative proton density of the macromolecular pool in a patient 
suffering from multiple sclerosis (MS). Figs. 5(a) and 5(b) are 
the images corresponding to the two acquisitions performed with 
and without the inversion pulse <x 3/ respectively. Fig. 5(c) shows 
the macromolecular fraction computed pixel -by-pixel according to 
Equation [17] . In brain tissue the macromolecular pool is asso- 
ciated to the myelin lipids and proteins and is therefore ex- 
pected to scale with myelin density. It can be nicely appreci- 
ated that white matter exhibits a higher macromolecular proton 
density than the gray matter does. MS plaques show a reduction 
in the macromolecular proton density. Of course, the macromo- 
lecular proton density in other tissues such as myocardial tis- 
sue or cartilage can be similarly determined using the tech- 
niques of the invention. 

Although exemplary embodiments of the invention have been 
described in detail above, those skilled in the art will appre- 
ciate that many additional modifications are possible in the ex- 
emplary embodiment without materially departing from the novel 
teachings and advantages of the invention. The pulse sequence 
illustrated in Fig. 4 represents only some of many possibilities 
for measuring the macromolecular proton density from the dilu- 
tion of phase modulated spins. In this sequence, a RF pulse (oc 3 ) 
is applied to disturb the equilibrium of the labelled spins by 
manipulating the spins of the free proton pool. Alternatively 
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disturbance of the equilibrium can be achieved with a sequences 
like in Fig. 4, but which uses single or trains with flip angles 
a 31 ....a 3n of off-resonant RF pulses instead of the third RF 
pulse of flip angle ot 3 in order to saturate the spins in the re- 
stricted proton pool. Such variation of the preferred pulse se- 
quence is illustrated in Fig. 6. Also, this principle can like- 
wise be applied to sequences like in Fig. 4, but where the flip 
angle of the inversion pulse deviates from 180° , or to a scheme 
were x 2 or x 3 are varied over successive scans. Accordingly, all 
such modifications are intended to be included within the scope 
of this invention as defined in the following claims. 



